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Abstract 

We discuss how the curvature and the strain density of the atomic 
lattice generate the quantization of graphene sheets as well as the 
dynamics of geometric quasiparticles propagating along the constant 
curvature/strain levels. The internal kinetic momentum of Rieman- 
nian oriented surface (a vector field preserving the Gaussian curvature 
and the area) is determined. 



1 Introduction 

Graphene and other "one-atom thick" giant 2D-molecules materialize in a 
sense the mathematical notion of abstract surfaces pQ. However, to be geo- 
metrically stable, materialized surfaces of this kind freely suspended in 3D 
Euclidean space, are to be curved and strained [2J-0- The curvature and 
strain generate some intrinsic fields which act on charge carriers similarly 
to magnetic and electric fields [8] [12]. The pseudo-magnetic field arising in 
this way forces the trajectories of charge carriers to form cycles analogous 
to the Larmor ones. These cycle currents can be considered as intrinsically 
generated "geometric" quasiparticles, whose size (de Broglie wavelength) cor- 
relates with the curvature radius and effective length of the lattice strain. 
Certainly, the size of geometric quasiparticles is smaller than the size of the 
geometrically stable area in which they live. 

Quasiparticles of this kind are able to propagate as a whole along the 
surface in the absence of any external fields, just due to the inhomogeneity 
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of the Riemann metric and strain. This dynamics, in the principle semiclas- 
sical approximation, preserves the state density as well as a curvature/strain 
symplectic form (Poisson brackets) on the surface [13]- [16]. At the quantum 
level, this leads to the appearance of a quantum structure making the surface 
coordinates to be noncommutative (like in the case of Landau-Peierls guiding 
center coordinates on the plane [H]-[IH]). The surface area has to be treated 
as a quantum "phase space", where the role of "Plank scale" is played by 
the inverse scale of the curvature or/and of the lattice strain density. Thus, 
ID -materials attaining geometric stability become quantum surfaces. 

Note that in general the strength of the intrinsic pseudo-magnetic field 
in graphene is composed of two sources as 



where the sing ± reflects the direction of the pseudospin, K is the Gaussian 
curvature of the graphene surface, and the function s can be called the strain 
density of the atomic lattice on the surface. The strain density is deter- 
mined by using a linear combination of first derivatives of the strain tensor 
components in a specific coordinate system attached to the lattice axes. 

The geometric quasiparticles on the graphene surface exist if and only 
if at least one of two magnitudes (1.1) does not vanish. The bands on the 
surface on which both the strain density s and the curvature K are zero or 
small are areas of quantum instability where the carbon 2D-lattice is flat, 
not stretched and therefore is not going to keep its surface geometry, but 
transforming to some different shape (tube, fullerene, schwarzite) or just 
crumpling to somewhat not two-dimensional. These unstable bands can be 
treated as "articulations" joining stable and quantized pieces of the graphene 
surface. 

From this viewpoint we considered in [TH] the charge carriers spectrum on 
graphene-like surfaces taking into account the strength of the external mag- 
netic field or the internal strain and ignoring the curvature contribution. In 
the present note, we complete this consideration by including the graphene 
curvature. We especially look at regions, where the curvature contribution 
dominates over the strain, and obtain the dynamics of geometric quasiparti- 
cles preserving the curvature and the area of the surface. The generator of 
this flow is a vector field which can physically be treated as an internal kinetic 
momentum of the Riemannian surface due to its curvature inhomogeneity. 

This classical picture of the geometric quasiparticle dynamics in graphene 
sheets is essentially corrected by the quantum topological condition a la 



s ± K/2, 
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Planck. The number of quantum states of the quasiparticle turns out to 
be proportional to the integral of density (1.1) over the grapheme area in 
question. In the case of small strain s ~ 0, as we demonstrate below by 
the Gauss-Bonnet theorem, the only way to have large enough number of 
quantum states, more than 1, is to assume that K < 0. Thus one can con- 
clude that not strained graphene areas of positive Gaussian curvature repel 
geometrical quasiparticles; these objects can naturally live in areas of negative 
curvature or, alternatively, they need a strong enough strain of the atomic 
lattice. 

2 Graphene algebra 

Charge carriers in graphene at energies near the bottom of the conductivity 
zone mimic the Dirac fermions [20] , [21] . The simplest version of the quantum 
Hamiltonian is the following (for details and generalizations, see, e.g., [6], |22j, 
121): 

H = vy ■ p, f~ 10 s cm/sec, (2.1) 

where 7 us a pseudospin and p is the kinetic momentum. In each local coor- 
dinate system q = (q 1 , q 2 ) on the graphene (orientable) surface, the following 
relations hold between components of 7 and p: 

[ 7 J ,7 m ] + = 2g jm (q), fe,7 m ] = itaybW, (2.2) 

and also relations involving coordinates: 

y,p m ] = mi, [q\q m ] = 0, ^, 7 m ] = 0. (2.3) 

Here [•, •] denotes the usual commutator, and [•,•]+ denotes the anticommu- 
tator. 

The mutual relations between components of the kinetic momentum p 
are the following |16j : 

\pj,p m ) = th 2 (S jm (q) + ^R s i jm {q)j sl ), (2.4) 

where 

l d = \b%l\ (2.5) 
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The tensor g^ 1 = {(g^ m )) in (2.2) represents the inverse metric on the surface. 
The tensor Rj ms i = gjrR r ms i in (2.4) is the curvature of the metric connection 
with the Christoffel symbols TJj: from (2.2): 



One can also include the torsion terms appearing due to dislocations in 
the atomic lattice [21] into connection coefficients T^. 

The tensor S = ([Sj m ]) in (2.4) is generated by the internal strain density 
s of the atomic lattice on the surface. More precisely, S has the form S = 
sJ, where the skew-symmetric covariantly constant tensor J is defined by 
the relation Ji 2 = ±\/det g, where the sign (±) detects the consistency or 
inconsistency of the given local coordinate system with the chosen orientation 
of the surface. 

Note that the Hamiltonian (2.1), in order to be self-adjoint has to be 
considered in the Hilbert space of half-densities over the given surface. This 
follows from the second relation in (2.2) and from the identities di In y^detg = 
Tjj, which are consequences of (2.7). 

Relations (2.2)-(2.4) generate an associate algebra indeed, since the Ja- 
cobi identities for double commutators and anticommutators hold. Namely, 
the Jacobi identity for the triple p,p,j follows from the definition (2.6) and 
the property R s ij m = Rjmsh the Jacobi identity for the triple p, 7, 7 follows 
because the connection (2.7) preserves the metric g = ([gj m ])] the Jacobi 
identity for the triple p, p, p follows from the closedness of the form 



and from the second Bianci identity (actually, in our two-dimensional case it 
holds automatically). 

Note that the commutators r y sl (2.5) occurring in (2.4), together with 7 m , 
generate the "Lie algebra" 



'msl ^s 1 lm ^l 1 - sm ' sk Im Ik sm) 



(2.6) 



(2.7) 



& = ( y l/2)S jm (q)dq j A dq m = ±s v / detgdq 1 A dq 2 



(2.8) 




2*arv-<rV), 

2i( g lm -f sl + g ms i jl + g S3 i lm + g jl i ms ). 



(2.9) 



Also note that 



(I s1 ) 2 



1 



(2.10) 



det g 
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This follows from the anticommutation relations (2.2). 

Our analysis of the algebra (2.2)-(2.4) is based, first of all, on the key 
relation (2.4). This shows that 

h 

p~- (2.11) 

t * 

where U is the characteristic scale of the strain-curvature field S+^Rj. Then 

W,p')~j, (2.12) 

by the first relation (2.3), where I stands for the characteristic scale of "large" 
inhomogeneity of the graphene sheet and q' = q/l and p' = (l*/h)p are 
normalized coordinates and momenta, q' ~ 1 and p' ~ 1. If the parameter 
e = l*/l in (2.12) is small, then one can separate the "slow variables" q' 
(for which [q',p'] ~ e) from the "fast variables" 7, p' (for which [7,7]+ ~ 1, 
\p',p'] ~ 1) by using the standard adiabatic approximation. 

The spectrum of the Hamiltonian (2.1) H = (Hv/l*)j ■ p' can be read- 
ily computed in the subalgebra of fast variables 7, p' producing a series of 
"Landau levels", i.e., energies of different size circular currents. Each "Lan- 
dau level" , except for the zero one, is actually a function depending on the 
slow variables. These variables are not just q', but are chosen from the ad- 
ditional condition that they commute with the fast variables up to e 2 . This 
condition can be achieved only if the slow variables are admitted to be non- 
commutative (the way in which the "leading center" coordinates appeared). 
Finally, one obtains a series of Hamiltonians over the surface with nontrivial 
commutators between coordinates. These Hamiltonians determine quantum 
states and the classical dynamics of geometric quasiparticles on the graphene 
quantum surface, which could be very useful in the "strain electronics" [25J. 

3 Strain quasiparticles 

Recall results of [16] for the case in which the strain dominates the curvature 
on the right-hand side of (2.4). 

Theorem 3.1. Assume that the tensor S = ([Sj m ]) does not degenerate on 
the area in question of the graphene surface and dominates the curvature 
field. Then the following statements hold. 
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(i) The Hamiltonian (2.1) in the low-energy approximation is equivalent 
to the direct sum of the "Landau level" Hamiltonians 

U{Q) = ±hvVk-47rN(Q), k = 0,1,2,... (3.1) 

where 

N = f (l/(2n)) v / \detS\/detg = \s\/2n (3.2) 

on the quantum surface with nontrivial commutation relations between coor- 
dinates: 

[Q j , Q m ] = iS~ ljm {Q) + small corrections. (3.3) 

The "small correction" summands in (3.3) are chosen to provide the correct 
behavior (invariance) after change of local quantum coordinates Qi in higher 
orders with respect to the semiclassical small parameter I* /I, where /* char- 
acterizes the scale of the strain density 

\s\ ~ I//* 2 . (3.4) 

(ii) For k ^ the classical dynamics of quasiparticles on the surface 
generated by the kth Hamiltonian (3.1) reads 

f = {n,Q} (3.5) 

with respect to the Poisson brackets {•, •} on the surface corresponding to 
relations (3.3), i.e., {A, B} = S^ m d m AdjB . The Hamilton-type system 
(3.5) reads as an equation of Maxwell-Lorentz type: 

4:71 

curlS* = ± — j, (3.6) 
where the "current density" j is defined by 

de f I 71 T^r^\3/2 d Q 



j 3 tfN(QY>^. (3.7) 

(iii) The function N (3.2) and the form & (2.8) are preserved by the flow 
generated by the dynamical system (3.6). The Planck-type discretization rule 
for the symplectic form (2.8) or, equivalently, the discretization rule for the 
integral strain 

1 f 1 

±— / sda = n+- iV = jV„, n = 0,1,2,..., (3.8) 



27r./E[jv] - 
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where dT\N] — {N(Q) = N} and da is the Riemannian meassure on the 
surface, implies the semiclassical asymptotics of the near-zero eigenvalues of 
the Hamiltonian (2.1): 

E Kn sa ±hvy/k ■ 4nN n . (3.9) 

Note that in view of (3.8) the function N (3.2) determines the state den- 
sitjQ of quasiparticles on "Landau levels" . Quasiparticles propagate along 
the curves {N(Q) = N n } surrounding the ares with the discrete flux (3.8) on 
the graphene surface. If the energy (3.9) is about the Fermi energy e f , the 
bound (3.4) for the strain density s correlates with the value of the Fermi mo- 
mentum pf = Sf/v and with the estimate for the kinetic momentum (2.11); 
the spatial quasiparticle size I* then correlates with the Fermi wavelength 
l F = h/p F . 



4 Curvature quasiparticles 

In contrast to the previous section, assume that the curvature contribution 
dominates the lattice strain contribution in (2.4). In such a situation, one 
can replace the commutation relation (2.4) of the graphene algebra by the 
relation 

[p i; p m ] = {ih 2 /A)R sljml sl . (2.4 a) 

The entire algebra (2.2), (2.3), (2.4a) is very interesting from the mathemat- 
ical point of view, because it is generated by the metric tensor exclusively. 
Therefore, any consequences derivable from representation theory and spec- 
tral theory for this algebra and for the Hamiltonian (2.1) contain information 
on geometrical properties of Riemannian surfaces. For instance, the proce- 
dure of adiabatic separation of variables and reduction to the "Landau lev- 
els", briefly described at the end of Sec. 2, produces a curvature preserving 
flow on the Riemannian surface which can physically be interpreted as the 
Hamiltonian flow of quasiparticles in graphene. Let us now go into details. 

The principal difference of the purely "curvature case" from the purely 
"strain case" is the presence of generators 7 s/ on the right-hand side of (2.4a). 
But since we deal with a 2-dimensional surface, the only nonzero generator 

1 Actual density of states in graphene equals 4iV because of two possible values of a 
pseudospin and two possible choices of a valley (corners of the Brillouin zone) . 
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is 7 12 = — 7 21 , and (2.4a) reads 

ih 2 

[pi,Pa] = ^-^i2i27 12 - 

It follows from (2.10) that the spectrum of 7 12 consists of numbers ±1 / y/det g 
at each point of the surface. Thus, on the eigenspaces of 7 12 , one obtains a 
scalar right-hand side in the relation for the momentum components, 

[pi,p 2 ] = ±"2" ^1212/A/deTg = ±— K y/defcg, 

where K stands for the Gaussian curvature of the surface, K = f \Rsij m g S3 g lm ■ 
Now we can just apply the results of [16] claimed in the previous section 
(Sec. 3) by choosing the value |if there instead of the strain density s. For 
instance, the function (3.2) is 

I if I 

N= 1 -^ (4.1) 
4tt V ; 

in this case. Thus, we obtain the following "curvature copy" of Theorem 3.1. 

Theorem 4.1. Let the surface be oriented. Assume that the Gaussian cur- 
vature K does not vanish on the area in question of the surface. Then the 
Hamiltonian (2.1) over the metric generated algebra (2.2), (2.3), (2.4a), in the 
low-energy approximation, is equivalent to the direct sum of the Hamiltonians 



H = ±fivy/k\K\(Q), k = 0,1,2,,... . (4.2) 

In (4.2), the noncommutative coordinates Q = (Q 1 ,^ 2 ) on the quantum sur- 
face obey the relation 



[Q 3 



, Q m ] = T i (J-J- 1 ^ (Q) + small corrections. (4.3) 



The pair of signs =F on the right-hand side of (4.3) corresponds to the pair of 
opportunities to choose the direction of the pseudospin (i.e., some eigensub- 
space of the generator 7 12 ) . The notion of "small corrections" was explained 
in Theorem 3.2 (a). 

For k 0, the classical dynamics of quasiparticles on the surface gener- 
ated by the kth Hamiltonian (4.2) and by relations (4.3) is given by 

mJQ = ±kh(J- 1 d\n \K\)(Q), (4.4) 
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where the effective mass m* is determined by the relation H = m^v 2 . 

The flow generated by the dynamical system (4.4) preserves the Gaussian 
curvature K and the surface area da. The Planck-type discretization rule, 
corresponding to relations (4.3) or, equivalently, the discretization rule for 
the integral curvature 

±± f Kfo = n+\ N = N n , n = 0,l,2,..., (4.5) 

def 

where dT,[N] = {\K\ = 47riV} 7 implies the semiclassical asymptotics (3. 9) for 
the near-zero eigenvalues of the Hamiltonian (2.1). 

Note that, by (4.2), one can estimate the nonzero "Landau levels", 

where /* stands for the effective curvature radius, K ~ Comparing 
with Ep = Hv/If, we see that the Fermi wavelength lp ~ l*/Vk~ correlates 
with the curvature radius. We also see from (4.5) that the quasiparticles 
propagate along the curves surrounding areas with discrete values of the 
integral Gaussian curvature. This dynamics is controlled by the system (4.4) 
whose right-hand side kHJ^dln \K\ can be referred to as the internal kinetic 
momentum of the Riemannian surfaceion the kih. "Landau level"). 

The dynamics of strain quasiparticles due to lattice stretch was described 
in (3.6) by using the electromagnetic terminology, including "current density" 
and "Maxwell-Lorentz equation" . In contrast to this approach, we represent 
the dynamics of curvature quasiparticles in a "mechanical" form (4.4) by 
introducing the notion of effective mass m* = -yJk\K\. However, it should 
be noted that equation (4.4) is not of Newton or Einstein type and, as seems, 
has no direct analogs in mechanics or general relativity theory. The flow 
generated by the internal momentum of the surface is due to the spinor 
framing (2.2), (2.3), (2.4a) originated from the metric field on the surface. 

Remark 4.1. Note that by the Gauss-Bonnet theorem the integral curvature 
of any piece S of surface, on which K > 0, is estimated as 

— / Kda<l, (4.6) 
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and the equality in (4.7) is realized only for the closed sphere (fullerene) 
£ ~ § 2 . Thus, the discretization rule (4.6) either does not hold or holds 
only for one number n = 0. This means that on such pieces of the graphene 
surface there cannot exist quantum states at all or only one state exists (for 
each Landau level k = 1,2,...). We conclude that the areas on the graphene 
sheet with K > repels the states of curvature quasiparticles^. 

But on pieces with K < 0, the discretization rule (4.6) can hold for many 
values of n, and so, many quantum states of the quasiparticles can exist 
in these areas. Therefore, the purely curvature quasiparticles, as quantum 
objects, naturally live in negatively curved graphene areas. 

Except trivial saddle surface or one-sheet hyperboloids (e.g., "worm- 
holes" |57])one can mention as interesting cases "schwartzites" [25], and 
the carbon foam [30]. For this type of surfaces £ the total integral curvature 
(in the compact case) is given by the Gauss-Bonnet: 



l -J tfd<7=l-f|(£), 



Air 

where q stands for the topological genus. Some phenomena observed in topo- 
logically complicated graphene-type structures [31] probably can be related 
to the existence of the curvature quasiparticles currents on negatively curved 
surfaces. 

Remark 4.2. In general situation, the "Landau-level" Hamiltonians in graphene 
sheets look as 



H = ±hc\l2k\s±^K\, k = 1.2 



where the mixture of the strain density and the curvature controls the dy- 
namics and the spectrum of geometric quasiparticles. The integral of the 
mixed magnitude \s ± \K\, presented in the general discretization rule simi- 
lar to (3.8) or (4.5), due to a large contribution of strain, can take values in a 
wide interval, and therefore many quantum states of geometric quasiparticles 
can exist even on positively curved graphene pieces (see, for instance, in [32"]). 
On the other hand, on negatively curved graphene sheets the contribution of 
strain is small enough [33J, and so the curvature effects probably dominate 
indeed. 



2 This is in good correspondence with the statement of [26 that pentagon rings in the 
carbon lattice repel the charge density (recall that pentagons are the cause of positive 
curvature and heptagons — of negative curvature in graphene). 



10 



References 

[1] J. C. Meyer, A. K. Geim, M. I. Katsnelson, K. S. Novoselov, T. J. Booth, 
and S. Roth, "The structure of suspended graphene sheets," Nature 446, 
60-63 (2007); arXiv: |cond-mat/0701379 

[2] P. E. Lammert and V. H. Crespi, "Geometrical perturbation of graphene 
electronic structure," Phys. Rev. B61, 7308-7311 (2000). 

[3] D. R. Nelson, T. Piran, and S. Weinberg (eds.) Statistical Mechanics of 
Membranes and Surfaces (World Scientific, Singapore, 2004). 

[4] S. Morozov, K. Novoselov, F. Scheldin, L. A. Ponomarenko, D. Jiang, 
and A. Geim, "Strong Suppression of Weak Localization in Graphene," 
Phys. Rev. Lett. 97, p. 016801 (2006). 

[5] A. F. Morpurgo and F. Guinea, "Intervalley scattering, long-range disor- 
der, and effective time-reversal symmetry breaking in graphene," Phys. 
Rev. Lett. 97, p. 196804 (2006). 

[6] J. Gonzales, F. Guinea, and M. A. H. Vozmediano, "The electronic 
spectrum of fullerenes from the Dirac equation," Nucl. Phys. B 406 
(3), 771-794 (1993); arXiv: |cond-mat/9208004[ 

[7] E.-A. Kim and A. H. Castro Neto, "Graphene as an electronic 
membrane," Europhys. Lett. 84 (5), p. 57007 (2008); arXiv: 
|cond-mat/0702562| r2 

[8] H. Suzuura and T. Ando, "Phonons and electron-phonon scattering in 
carbon nanotubes", Phys. Rev. B65, 235412 (2002). 

[9] K. Sasaki, Y. Kawazoe, and R. Saito, "Local energy gap in deformed 
carbon nanotubes", Prog. Theor. Phys. 113, 463-480 (2005). 

[10] A. Cortijo and M. A. H. Vozmediano, "Electronic properties of curved 
graphene sheets", Europhys. Lett. 77, p. 47002 (2007); 
M. A. H. Vozmediano, F. de Juan, and A. Cortijo, "Gauge fields and 
curvature in graphene," J. Phys. Conf. Ser. 129, p. 012001 (8pp.) (2008); 
arXiv: cond-mat/0807.3909vl. 



11 



[11] F. Guinea, A. K. Geim, M. I. Katsnelson and K. S. Novoselov, "Gener- 
ating quantizing pseudomagnetic fields by bending graphene ribbons", 
Phys. Rev. B81, 035408 (2010); 

M. A. H. Vozmediano, M. I. Katsnelson, and F. Guinea, "Gauge fields 
in graphene" ; arXiv: cond-mat/1003.5179v2. 

[12] N.-C. Yeh, M.-L. Teague, S. Yeom, B. L. Standley, R. T.-P. Wu, 
D. A. Boyd, and M. W. Bockrath, "Strain- induced pseudo- magnetic 
fields and charging effects on CVD-grown graphene"; arXiv:cond-mat, 
1009.0081vl. 

[13] M. V. Karasev, "Magneto-metric Hamiltonians on quantum surfaces in 
the configuration space," Russ. J. Math. Phys. 14 (1), 57-65 (2007). 

[14] M. V. Karasev, "Integral geometric current, and the Maxwell equa- 
tions as a Hamiltonian system on configuration surfaces," Russ. J. Math. 
Phys. 14 (2), 134-141 (2007). 

[15] M. V. Karasev, "Geometric dynamics on quantum nano-surfaces and 
low-energy spectrum in a homogeneous magnetic field," Russ. J. Math. 
Phys. 14 (4), 440-447 (2007). 

[16] M. V. Karasev, "Quantum geometry of nano-space," Russ. J. Math. 
Phys. 15 (3), 417-420 (2008); "Magneto-metric current competing with 
Hall current on curved nanosurfaces," Proc. of Intern. Forum on Nan- 
otechnologies, Moscow, 2008, Vol. 1. p. 31 (in Russian). 

[17] R. Peierls, "Zur Theorie des Diamagnetismus von Leitungselektronen," 
Z. Phys. 80, 763-791 (1933). 

[18] R. Jackiw, "Physical instances of noncommutative coordinates," Lecture 
Notes in Phys. 616, 294-304 (2003); arXiv: |hep-th/0110'057l 

[19] G. Dunne and R. Jackiw, "Peierls substitution and Chern-Simons quan- 
tum mechanics," Nuclear Phys. B. Proc. Suppl. 33C, 114-118 (1993). 

[20] P. R. Wallace, "The band theory of graphite," Phys. Rev. 71, p. 622 
(1947). 

[21] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, 
I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, "Two-dimensional gas 
of massless Dirac fermions in graphene", Nature 438, 197-200 (2005). 



12 



[22] M. S. Foster and I. L. Aleiner, "Graphene via large N: A renormalization 
group study," Phys. Rev. B77, p. 195413 (2008). 

[23] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and 
A. K. Geim, "The electronic properties of graphene," Rev.Mod. Phys. 
81, 109-162 (2009); arXiv: cond-mat/0709.1163v2 

[24] F. de Juan, A. Cortijo, and M. A. H. Vozmediano "Dislocations and 
torsion in graphene and related systems," Nuclear Phys. B 828 (3), 
625-637 (2010); arXiv: cond-mat/0909.4068vl. 

[25] T. Low and F. Guinea, "Strain-induced pseudomagnetic field for novel 
graphene electronics", Nano Lett., 10 (9), 3551-3554 (2010). 

[26] S. Azevedo, C. Furtado, and F. Moraes, "Charge localization around 
disclintions in monolayer graphene", Phys. Status Solidi (b) 207 (2), 
387-391 (1998). 

[27] A. L. Mackay and H. Terrones, "Diamond from graphene" , Nature 352, 
762 (1991); 

E. R. Margine, A. N. Kolmogorov, D. Stoikovic, J. O. Sofo, and 
V. H. Crespi, "Theory of genus reduction in alkali-induced graphiti- 
zation of nanoporous carbon", Phys. Rev. B 76, 115436 (2007). 

[28] D. Vanderbilt and J. Tersoff, "Negative-curvature fullerene analog of 
C 60 ", Phys. Rev. Lett. 68 (4), 511-513 (1992). 

[29] H. Terrones, M. Terrones, E. Hernandez, N. Grobert, J.-C. Charlier, and 
P. M. Ajayan, "New metallic allotropes of planar and tubular carbon", 
Phys. Rev. Lett. 84, 1716-1719 (2000). 

[30] A. V. Rode, E.G. Gamaly, and B. Luther-Davies, "Formation of claster- 
assembled carbon nano-foam by high repetition-rate laser ablation", 
Appl. Phys. A70, 135 (2000). 

[31] V. I. Tsebro and O. E. Omel'yanovskii, "Persistent currents and cap- 
ture of magnetic flux in multiconnected carbon nanotubular structure" , 
Uspekhi Fiz. Nauk 170 (8), 906-912 (2000); 

N. Park, M. Yoon, S. Berber, J. Ihm, E. Osawa, and D. Tomanek, "Mag- 
netism of all-carbon nanostructures with negative Gaussian curvature" , 
Phys. Rev. Lett. 91 (23), 237-204 (2000). 



13 



[32] N. Levy, S. A. Burke, K. L. Meaker, M. Panlasigui, A. Zettl, F. Guinea, 
A. H. Castro Neto, and M. F. Crommie, "Strain-induced pseudo- 
magnetic fields greater than 300 Tesla in graphene nanobubbles" , Sci- 
ence 329 (5991), 544-547 (2010). 

[33] T. Lenosky, X. Gonze, M. Teter, and V. Elser, "Energetics of negatively 
curved graphitic carbon", Nature 355, 333 (1992). 



14 



